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Recently, a growing interest has been addressed to the electrical properties of bacteriorhodopsin 
(bR), a protein belonging to the transmembrane protein family. Several experiments pointed out 
the role of green light in enhancing the current flow in nanolayers of bR, thus confirming potential 
applications of this protein in the field of optoelectronics . By contrast, the mechanisms underlying 
the charge transfer and the associated photocurrent are still far from being understood at a mi¬ 
croscopic level. To take into account the structure-dependent nature of the current, in a previous 
set of papers we suggested a mechanism of sequential tunneling among neighbouring amino acids. 

As a matter of fact, it is well accepted that, when irradiated with green light, bR undergoes a 
conformational change at a molecular level. Thus, the role played by the protein tertiary-structure 
in modeling the charge transfer cannot be neglected. The aim of this paper is to go beyond previous 
models, in the framework of a new branch of electronics, we called proteotronics, which exploits the 
ability to use proteins as reliable, well understood materials, for the development of novel bioelec- 
tronic devices. In particular, the present approach assumes that the conformational change is not 
the unique transformation that the protein undergoes when irradiated by light. Instead, the light 
can also promote a free-energy increase of the protein state that, in turn, should modify its internal 
degree of connectivity, here described by the change in the value of an interaction radius associated 
with the physical interactions among amino acids. The implemented model enables us to achieve a 
better agreement between theory and experiments in the region of a low applied bias by preserving 
the level of agreement at high values of applied bias. Furthermore, results provide new insights on 
the mechanisms responsible for bR photoresponse. 


I. INTRODUCTION 

Bacteriorhodopsin (bR) is the best known protein in 
the family of opsins, proteins conjugated with a molecule 
of retinal and able to convert visible light into electro¬ 
static energy [l|. This protein is found in a primeval 
organism, the Halobacterium salinarum, specifically in a 
part of its cell membrane called purple membrane (PM), 
since its color. This membrane, 5 nm thick, is a natural 
thin film, essentially constituted by few lipids and these 
proteins organized in an hexagonal lattice 0. 

A large number of studies has been carried out on bR 
in the field of biophysics and physicochemistry , and 
many aspects have been unveiled. As relevant examples 
we cite: (i) the photoinduced isomerization of the retinal 
embedded in bR Q, (ii) the conformational change of 
bR associated with the retinal isomerization Q , (iii) the 
importance of environmental conditions in the photocycle 
development [1,0. 

Patches of PM have been used for several purposes: to 
produce metal-protein-met al j unctions |10l - ll2j| . to per¬ 
form c-AFM investigations [la [l4l| . to develop solar cells 


* eleonora.alfinito@unisalento.it, http://cmtgl.unile.it/eleonoral.html 
^ lino.reggiani@unisalento.it 


of new generation [Tl|, etc. As a matter of fact, films of 
bR resist to thermal, electrical and also mechanical stress 
[lol - [l^ and show a substantial photocurrent when irra¬ 
diated by a visible (green) light EOill. Therefore, 
bR can be used as an optoelectrical switch, to convert 
radiant energy into electrical energy [l^ . in pollutants 
remediation systems El: to produce optical memories 
[10, to control neuronal and tissue activity El, [20 , etc. 

The commonly accepted view concerning the protein 
activation is the following: A photon is absorbed by the 
retinal molecule contained in each protein, then caus¬ 
ing the bending of this molecule. As a consequence, the 
protein undergoes a change of its tertiary structure, fol¬ 
lowing a cycle of transformations that arrives to release a 
proton outside the cell membrane. Finally, reprotoniza- 
tion of the retinal molecule by Asp96 restores the native 
configuration. Some crystallographic investigations have 
been performed on this protein to determine its config¬ 
uration in the different steps of the cycle [0. This is a 
particularly hard task, since the X-ray radiation could 
modify the protein structure, and only recently the puz¬ 
zle of many contradictory results starts to be recomposed 
[10. At present, a rather complete description of the 
protein is given only for the native and the active L-state 

Measurements of the protein current-voltage (I-V) 
characteristics were reported in several papers [l0l - [l4l| . 
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To this purpose, samples made of patches of PM were 
anchored on a conductive substrate and connected to an 
external circuit. The connection was made with : i) an 
extended traiisparent conductive contact [iM3> ii) a tip 
of a c-AFM [il, [ 13 . In both cases, the measured cur¬ 
rent was found to be quasi-Ohmic at the lowest bias, and 
strongly super-linear at increasing bias. Furthermore, 
when the sample was irradiated with green light 
a significant photocurrent was observed. There was a 
clear proof that the charge transfer is mainly due to the 
protein [13, [13 ■ There was also a high resistance channel 
due to the lipid membrane, which is detected in exper¬ 
iments involving a membrane deformation fl^ . 1^ . In 
the absence of membrane deformation, this channel can 
be neglected. The charge transfer through the protein 
was attributed to a tunneling mechanism [^ . 

In particular, the presence of a current (in dark) well 
above possible leakage components, and of a photocur¬ 
rent (in light) supports the hypothesis of a mechanism 
of charge transfer intrinsically dependent on the protein 
tertiary structure 

Since a long time, the interaction of electromagnetic 
fields with biological matter is the object of many in¬ 
vestigations, mainly for the damages produced by ioniz¬ 
ing radiations. As far as known, sunlight that reaches 
the Earth is largely composed of non-ionizing radiations 
whose main effect on biological matter is heating. In par¬ 
ticular, for proteins, this should lead to a global energy 
enhancement, regardless of the protein specific conforma¬ 
tional state, as confirmed by recent experiments showing 
the critical role of temperature in current measurements 
[T3 - Therefore, we conjecture that in a sample of pro¬ 
teins, like a patch of purple membrane, light gives rise to 
different effects. From one side, the retinal modification 
with the consequent conformational change from the na¬ 
tive to the active state, from another side, a net transfer 
of energy to the whole protein with a consequent increase 
of its free-energy. As a general issue, both these effects 
should contribute to the protein activation. 

The present paper addresses this issue by account¬ 
ing simultaneously for these two effects in a computa¬ 
tional/theoretical model called INPA (impedance net¬ 
work protein analogue). This approach describes the 
electrical characteristics of a protein by using a network 
of impedances. In previous investigations, the local inter¬ 
action of a photon with the retinal has been investigated 
by considering the corresponding change of the network 
structure The novelty of the present paper consists in 
the further introduction of a global energy increase of 
the the whole protein due to the incident light. This is 
described by a change of the network connectivity both 
of the native and the active state. 

The methodological approach we follow points to the 
integration of different disciplines (molecular biology, 
physics, electronics) to develop a new generation of elec¬ 
tronic devices within a nano-bio-technology. This inter¬ 
disciplinary approach is leading to an entirely new disci¬ 
pline which we christen proteotronics ■ 


The paper is organized as follows. Section II summa¬ 
rizes the main steps of the INPA model and describes 
the improvements introduced on the basis of the dynam¬ 
ical evolution of the protein energy landscape. Section 
III reports and discusses the main results and suggests 
the opening of new perspectives. Major conclusions are 
summarized in Section IV. 


II. THEORETICAL MODEL 

The INPA model is based on a percolative approach 
that describes the protein like a network of links and 
nodes. A node represents a single amino acid and its spa¬ 
tial position is the same of the corresponding Ca atom. 
A link joins a couple of nodes, and represents the interac¬ 
tion between amino acids [334^ . The protein structure 
in its native or active state is taken by public databases or 
homology modeling [H, [s^ , thus the node configuration 
reproduces the protein backbone. Then, couples of nodes 
are connected with the rule that they must be be closer 
than an assigned interaction radius. Re- In this way, the 
number of links, N, depends on the value of Rc and is 
in the range 0 < Naa < Naa{Naa - l)/2, with Naa the 
number of amino acids pertaining to the given protein. 
In the present case, the macroscopic quantity of interest 
is the static I-V characteristic. Therefore, the network is 
drawn like an electrical circuit where an elementary re¬ 
sistance, is associated with each link between nodes 
i and j. Explicitly: 


where Aij = 7r(i?^ — lijjA), is the cross-sectional area 
between two spheres of radius Rc centered on the z-th 
and j-th node, respectively; hj is the distance between 
the sphere centers, p is the resistivity. 

By positioning the input and output electrical con¬ 
tacts, respectively, on the first and last node (more struc¬ 
tured contacts can be envisioned) for a given applied bias 
(current or voltage operation modes according to conve¬ 
nience) the network is solved within a linear Kirchhoff 
scheme and its global resistance, R, is calculated [23 - 
[13 • Accordingly, this network produces a parameter- 
dependent static I-V characteristic for the single protein, 
based on the standard relation: 

V = RI. (2) 

To account for the super-linear behaviour of current at in¬ 
creasing voltages, a tunneling mechanism of charge trans¬ 
fer is included. In doing so, a stochastic approach within 
a Monte Carlo scheme [Ij, [13 [13, [234^ i s used. In 
particular, following the Simmons model |38l |. a mecha¬ 
nism containing two possible tunneling processes, a direct 
tunneling (DT) at low bias, and a Fowler-Nordheim tun¬ 
neling (EN) at high bias, is introduced. Therefore, the 
resistivity value of each link is chosen between a low value 
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Pmin, taken to fit the current at the highest voltages, and 
a high value p{V), which depends on the voltage drop 
between network nodes as: 


piV) = pmax (eV < $), (3) 

P(V) = PMAX(^) + Pmin(l “ (eV > $) (4) 

where pmax is the maximal resistivity value taken to 
fit the I-V characteristic at the lowest voltages (Ohmic 
response) and <i> is the height of the tunneling barrier be¬ 
tween nodes. The transmission probability of each tun¬ 
neling process is given by [H, [ 23 : 


pDT 


exp 





pFN 


= exp 


—a 





( 5 ) 

( 6 ) 


where Vij is the potential drop between the couple of 

i,j amino acids, a = and m is the electron 

effective mass, here taken the same of the bare value. 
The DT superscript refers to the low-bias, quasi-Ohmic 
response and the FN subscript refers to the high-bias, 
super-Ohmic response. 


III. RESULTS AND DISCUSSION 

By construction, both the current response at very low 
and very high bias exhibit an Ohmic behaviour with val¬ 
ues of the corresponding resistance differing for several 
orders of magnitude. This model was successfully used 
to reproduce the experiments of Ref. [l^ . The inputs pa¬ 
rameters were i?c=6 A, $=219 meV, pmax = 4 x 10^^ O 
A, for the low field resistivity and Pmin = 4 x 10"* OA for 
the high field resistivity. The protein tertiary structure 
was taken from the protein database [1^ , specifically the 
2NTU entry, an X-ray crystallographic measurement for 
bR native-state. 

The agreement between calculations and experiments 
was found to be satisfactory, also reproducing the cur¬ 
rent modifications due to the membrane indentation by 
the c-AFM tip [ 2 ^. On this basis, we found reasonable 
to take the same input to fit the response of the protein 
in light. At present, the only crystallographic entry de¬ 
scribing the complete protein in an active state is 2NTW, 
which gives account of the L-state (henceforth called the 
active state) of bR. This state is sensitive to the 550 nm 
light and precedes the M state (410 nm), which corre¬ 
sponds to a proton releasing. In the TV measurements, 
the proton releasing was not monitored and the current 
measured was only attributed to electron transfer M- 
When the activated configuration was used as input to 
fit the current response in the presence of light, the agree¬ 
ment with experiments was less satisfactory than that in 
dark. A possible way to overcome this drawback is to 
assume that the presence of light modifies not only the 



FIG. 1. Relative resistance variations vs the interaction ra¬ 
dius, Rc, for bR in the native and active state. The ellipse 
indicates the region of Rc values whose trend is in agreement 
with experiments. 


protein structure but also its connectivity properties. In 
the INPA model this modification is accounted for by 
changing the value of the interaction radius. To this pur¬ 
pose, Fig. [T] reports the role of the interaction radius in 
the calculation of the resistances of the native and active 
states. Numerical data are obtained at very low volt¬ 
ages, where the Ohmic regime strictly holds, and reflect 
the protein topology. The main results of Fig. [T] are: 

(i) the general low resolution between these states and, 

(ii) the presence of two regions in which the resolution is 
best appreciable, around Rc = 6 A and Rc = 9 A. The 
experiments are in agreement with Rc = 6 A. 

In the following the protein current responses are nu¬ 
merically analyzed for several values of Rc around 6 A. 
In particular, simulations for the single protein are per¬ 
formed for the three values, Rc= 5.8, 6.0 and 6.3 A and 
for both native and active state. Results are reported 
in Fig. [5] with panel (a) reporting the experimental data 
carried out in a bR macroscopic sample |10l | in dark and 
light. For a given protein state. Fig. shows a current 
enhancement by increasing Rc from 5.8 to 6.3 A. Fur¬ 
thermore, at increasing Rc, the differences between the 
current response in dark and light are more and more 
marked. 

The above results suggest that the activation mech¬ 
anism of a macroscopic sample of bR can be described 
within the single protein model by using: i) a confor¬ 
mational change (from 2NTU to 2NTW structure), ii) 
a connectivity change (i.e. a variation of the network 
interaction radius). 

More specihcally, we can envision a twofold mechanism 
of photon absorption: by the retinal, and by the whole 
protein. The former is responsible of the conformational 
change, the latter of a global energy increase of the pro¬ 
tein. Notice that, according with experiments [l^, the 
global energy increase is coherently used by the protein 
sample in enhancing its photocurrent response. In other 
words, the electromagnetic radiation impinging on the 
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FIG. 2. I-V characteristics of bR. Panel (a) reports the exper¬ 
iments carried out on nanolayers samples [1^ Panels (b,c,d) 
refer to data calculated within the INPA model (single pro¬ 
tein) for the different values of Rc reported in the figures. 
Symbols refer to numerical calculations, lines are guides to 
the eyes. Dashed lines and superimposed symbols refer to the 
active state; continuous lines and superimposed symbols to 
the native state. For 7?c=5.8 A the I-V characteristics are 
found to coincide for native and active states. 



FIG. 3. Schematic representation of the energy evolution by 
absorption of photons in a bR molecule. Upper panel depicts 
the conformational change from the native state U to an ac¬ 
tive state W induced by the absorption of a photon by the 
retinal and a successive global energy increase from W to W' 
induced by the absorption (from the whole protein) of other 
photons . Lower panel depicts an alternative possibility when 
the absorption of photons induces a global energy increase 
from U to U' of the native state and a successive absorption 
process induce a conformational change from U' to the energy 
level W' of the active state. 



FIG. 4. Schematic representation of the free-energy funnel 
landscape for native and active states. The minima of the 
free-energy corresponds to a minimum of the conformational 
entropy that measures the number of available microscopic 
states. The increase of connectivity at increasing free-energy 
is depicted by the increase of the interaction radius R^. 



FIG. 5. Gurrent fluctuations of bR in its active state at dif¬ 
ferent values of Rc 


protein may anyway produce the global effect of an en¬ 
ergy gain, while the local interaction of a photon with 
the retinal triggers the conformational change when the 
protein is in its native state. These mechanisms associ¬ 
ated with photon absorption are schematically depicted 
in Fig. [3l From one hand, when the native state, say 
U, becomes an active state, say W, a further irradiation 
should enhance the global energy of the active state. In 
this way, the active state is promoted to an upper energy 
value W. From another hand, when the U state does 
not undergo a conformational change, anyway its energy 
level can be promoted to an upper value U'; a further 
dose of light may drive this state to an active state W. 

Among the different ways used to describe the protein 
energy landscape at different stages of the folding one of 
the most accepted is the rugged funnel-diagram [^. In 
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this diagram, the protein folds from the molten state to 
the native (stable) state following many possible folding 
routes toward the minimum of a funnel-like energy sur¬ 
face. When the protein runs down in the energy funnel, 
it loses the spurious bonds and enforces those stabiliz¬ 
ing the minimal-energy configuration [s^. In doing so, 
it also reaches the minimum of the configurational en¬ 
tropy. Furthermore, the phase transition from a stable 
state at low energy to a stable state at higher energy is 
depicted in terms of a tunneling between the minima of 
a multiwell energy landscape [4^. As energy increases, 
the spurious connections do again appear and the protein 
can explore more microstates. In a very schematic way, 
this mechanism is pictured in Fig. |4] where a couple of 
funnels representing the native and an active state are 
superimposed. The conformational change corresponds 
to the transition from a funnel to another one. The min¬ 
imal energy between the two funnel stable minima is of 
the order of the eV. Otherwise, by rising the energy of 
the protein in the native state, it is possible to reach the 
overlapping region of the two funnels. Here the transi¬ 
tion from the native to the active state can occur without 
energy supply. 

Within the INPA model, the mechanism of energy in¬ 
crease is described by an increase of the Rc value. As a 
consequence, the network becomes more connected which 
implies an increase of the pathways for tunneling and of 
the number of possible current channels. This, in turn, 
leads to an increase of the instantaneous current fluc¬ 
tuations, as reported in Fig. [S] Here, the current fluc¬ 
tuations observed from simulations are reported for the 
active state, with an applied bias of 0.75 V and for two 
Rc values of 5.8 and 6.3 A. 

Following this scheme, the current response of samples 
made by monolayers of bR has been fitted by using a 
binary mixture of native and active states; the percent¬ 
ages of each state being a function of the Rc valu e 1^ . 
Specifically, a good fit of the experimental data [10l| is 
obtained by using: (i) for the sample in light, Rc = 6.3 
A and a binary mixture of 96% of 2NTW and a 4% of 
2NTU; (ii) For the sample in dark Rc = 5.8 A and 100% 
of 2NTU. 

The c-AFM experiment, performed in the absence of 
direct light [l^ , was previously fitted within a very good 
accuracy on the full bias range by using Rc = 6.0 A and 
100% of the 2NTW native state [2^. Since in these ex¬ 
periments one cannot exclude the presence of a certain 
amount of proteins IN the active state, in agreement with 
a value of Rc larger than the threshold value Rc = 5.8 
A here the fit with experiments is tested by using binary 
mixtures with an increasing percentage of active states. 
The fit is found to be sufficiently accurate with a per¬ 
centage of active state not larger than 40%. Figure [5] 
reports: (i) the experimental data [iSlil and, (ii) the 
single protein data rescaled by using the formula 

~ A (/n * '^n T la * ria) (7) 

where Is indicates the sample current, A is a numerical 



FIG. 6. I-V characteristics calculated by using different mix¬ 
tures of bR native and active states. Full circles refer to the 
mixture 4% of native states and 96% of active states with 
Rc=6.3 A. Open squares refer to a pure native state with 
7?c=5.8 A. Continuous (dashed) line refers to experimental 
data m in dark (light) in the bias range 0.1 1 V. In the 

inset the continuous line refers to experimental data in dark 
0. in the bias range 0.1 -;- 10 V. Circles refer to data calcu¬ 
lated with the pure native state with Rc=6.0 A, squares refer 
to data calculated with the mixture of 60% of native states 
and 40% of active states, with Rc=6.0 A. 


constant of the order of 10^ used to scale the single pro¬ 
tein current to the macroscopic data, In/a is the current 
of the single protein calculated with the native/active 
configuration, n„/Q is the fraction of native/active pro¬ 
tein expected in the sample. Of course, for a pure state, 
this formula reduces to the simple proportional rescaling: 

Is = Aln/a ( 8 ) 


Figure [7] reports the concentration of active 2NTW 
states in the samples vs the corresponding Rc values to be 
used in simulations. Symbols refer to values used in sim¬ 
ulations and the dashed curve is a fitting obtained from 
a sigmoidal Hill-like function that is commonly used in 
biochemistry to describe the percentage of proteins acti¬ 
vated by a ligand. Its validity in fitting several different 
physicochemical reactions is well known (4ll| . and writes: 


fix) 


x°‘ 

b + x°‘ 


X 100, 


X = 


Rc — Ro 
Rq 


(9) 


where / is the percentage of proteins in the active state, 
and i?o=5.8 A. For = b, half of the proteins in the 
sample have changed their configuration. Here, the best 
fitting parameters are a = 3.91 and b = 2.87 x 10“®, 
i.e. = (6.02 — i?o) /Ro- The full circles reproduce 
the experiments reported in Fig. [T] For the case of the 
experiments in Ref. [l^ . further binary mixtures with 
the percentages suggested by Eq. [S] (open circles) have 
been tested to to be consistent with experiments but to 
a less quantitative resolution of the photocurrent. 
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FIG. 7. Per cent concentration of proteins in the active state 
for the different values of the interaction radius to be consid¬ 
ered in simulations. Full circles refer to the values used in 
Fig. [6l dashed line is the fitting obtained with Eq. 0 empty 
circles refer to values intermediate between the experiments 
reported in Ref, m , see text. 


In the present context, the meaning of the function 
/ is the following: for an increasing number of photons 
impinging onto the sample, the free energy of the protein 
grows and, as a consequence, the value of Rc also grows. 
With Rc also the percentage of proteins moved to the 
active state grows because some of the photons hit the 
retinal. Finally, for Rc larger than about 7 A all the 
proteins in the sample are in the active state. Further 
amount of photons may only improve the free energy of 
the protein in the active state and the internal degree of 
connections. 


IV. CONCLUSIONS 

The paper investigates the mechanisms responsible for 
the photocurrent exhibited by monolayer samples of bac- 
teriorhodopsin in the presence of an impinging green 
light. To this purpose, use is made of the INPA model 
implemented to account for the change of connectivity 
of the single protein associated with the presence of the 
light. Previous results provided a satisfactory interpreta¬ 
tion of a set of accurate measurements, performed with 
the c-AFM technique, in nanolayer samples and in the 
absence of direct light. Accordingly, experiments were 
interpreted on the basis of the tertiary structure associ¬ 
ated with the native state of the single protein. However, 
a less satisfactory agreement was obtained in the region 
of low voltages when the same approach was applied to 


the case of monolayer samples of bR in the presence of 
light, and thus taking the tertiary structure of the pro¬ 
tein in its active state. To overcome this drawback, here 
we consider also the change in the connectivity of the 
protein state consequent to the enhancement of the free- 
energy level of the single protein induced by the presence 
of the light. The increase in connectivity is accounted for 
by an increase of the value of the interaction radius, Rc, 
already introduced to correlate the electrical properties 
with the tertiary structure of the protein. Accordingly, 
the new model interprets the photocurrent using a binary 
mixture of results pertaining to the native and active 
structures with the proper values of Rc [l^. Specifically, 
a satisfactory fit of experiments on nanolayers is ob¬ 
tained by using: (i) for the sample in light, Rc = 6.3 
A and a binary mixture of 96% of 2NTW and of 4 % 
of 2NTU, (ii) for the sample in dark, Rc = 5.8 A and 
100 % of 2NTU. The c-AFM experiments, performed in 
the absence of direct light, are quite finely reproduced by 
using a binary mixture containing up to 40% of 2NTW 
and i?c=6.0 A (see Fig. O. Therefore, the implemented 
model enables us to achieve a better agreement between 
theory and experiments in the region of low applied bias 
and does not modify previous findings at high values of 
applied bias. 

We notice that the process of protein activation, in 
particular for opsins, is still a very open topic [Ulil 
and the present approach aims to provide a further step 
for a better understanding of the subject. Environmen¬ 
tal effects, different from the presence of light, like tem¬ 
perature, the value of the pH, etc, should be responsi¬ 
ble for other activation mechanisms. Accordingly, more 
experiments, and structural information are necessary, 
and present results should give a further motivation to 
stimulate new experiments and formulate new theories. 
Finally, this research exploits the trend in which differ¬ 
ent emerging disciplines can converge in a new branch 
of science, we recently introduced as proteotronics [s^- 
Indeed, proteotronics aims to develop new devices based 
on the sensing properties of proteins. In doing so, protein 
responses to external stimuli have chances to be better 
understood and used to devise biodevices of relevant im¬ 
portance in applied sciences. 
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